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ABSTRACT 



Aims. Our goal is to find a minimalistic model describing the luminosity function and large-scale clustering bias of the 
X-ray selected AGN in the general framework of the concordance ACDM model. 

Methods. We assume that a simple population-averaged scaling relation between the AGN X-ray luminosity Lx and 
the host dark matter halo mass Mh exists. With such a relation, the AGN X-ray luminosity function can be computed 
from the halo mass function. Using for the latter the concordance ACDM halo mass function, we obtain the Mh — Lx. 
relation required to match the redshift dependent AGN X-ray luminosity function known from X-ray observations. 
Results. We find that with a simple power law scaling A/h oc ^x i '^^^ model can successfully reproduce the observed X- 
ray luminosity function. Furthermore, we automatically obtain predictions for the large-scale AGN clustering amplitudes 
and their dependence on the luminosity and redshift, which seem to be compatible with AGN clustering measurements. 
Our model also includes the redshift-dependent AGN duty cycle which peaks at the redshift z ~ 1 and its peak value 
is consistent with unity, suggesting that on average there is no more than one AGN per dark matter halo. For a typical 
X-ray selected AGN at z ~ 1, our best-fit Mh — Lx scaling implies low Eddington ratio I/x/i/Edd ~ 10~* — 10"'' and 
correspondingly large mass growth e-folding times, suggesting that the typical X-ray AGN are dominantly fueled via 
relatively inefficient 'hot-halo' accretion mode. 
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1. Introduction 

A realization that every galaxy of significant size har- 
bors a supermassive black hole (SMBH) (with masses ^ 
10^ — 10^ Mq) is probably one of the most remarkable dis- 
coveries of modern astrophysics. The growth of SMBHs, as 
manifested by the active galactic nuclei (AGN), has been 
observed over a broad range of electromagnetic energies, 
from radio to hard X-rays and gamma-rays. The observed 
AGN X-ray emission, believed to originate from the upscat- 
tering of softer accretion disk photons via inverse Compton 
mechanism by the hot electron corona, has proven to be the 
most effective way of selecting large samples of AGN over 
large cosmological volumes [^ This efficiency owes to the 
fact that in order to absorb X-rays one needs significant 
column densities of absorbing material, which is particu- 
larly true in the hard X-ray band, where iVn ^ 10^^ cm~^ 
is required. 

With XMM-Newton and Chandra X-ray observato- 
ries, two of the most advanced X-ray instruments in exis- 
tence, more than 20 deep extragalactic X-ray surveys over 
varying sky areas and limiting sensitivities have been per- 

^ In comparison, the deepest optical spectroscopic surveys 

typically give a factor of ~ 10 times less AGN per deg~^, and 

only utradeep optical variability studies arc able to generate 

comparable AGN s ky densities (e.g.,|Brandt & Hasinger 20051). 

http : //xmm . e sac.esa. int _ 



formed (Brandt & Hasinger 2005). The most noticeable 
amongst th ese are ^ 2 Ms Chandra Deep Field North 
([Alexander et al. 20031) over ~ 448 arcmin^ sk y area, ~ 4 
Ms Chandra Deep Field South (iXue et al.|2011|) over ~ 465 



and ^ 3 Ms XMM -Newton deep survey covering 
Chandra Deep Field South ( Comastri et al.|[2011 ). 

Despite the small sky areas covered by those deep sur- 
veys, a relatively large number of detected AGN has al- 
lowed reasonably good determination of the AGN lumi- 
nosity function (LF) and its evolution over cosmic tim e 



(Ueda et al. 2003 Hasinger et al. 2005 Aird et al. 2010) 



Even though these narrow survey geometries are not very 
suitable for measuring the spatial clustering properties 
of AGN, somewhat wider and shallower surveys like X- 
Bootes, XMM-LSS, AEGIS, and XMM-COSMOS have al- 
lo wed measure ment of th e two-po int clu stering statistics 



Coil et al. 2009 



(M urray eta l. 2005; Gandhi ct al. 2006 
GiUi et al.|2009; |Allevato et al.,2011^ , al'B^e'it with relatively 
large uncertainty on large scales. The large-scale clustering 
measurements are expected to improve significantly with 
the upc oming Spectrum- X-G amma/eROSITA n [j space 
mission (Predehl et al. 20101, which is planned to cover 
entire sky down to the limiting se nsitivity of ^ 10"^"' 
erg/s/cm^ in the 0.5 — 2 keV band (Merloni et al 
Kolodziget aLl[2012l ). 



2012 



ihttp : //chandra . harvard . edu 



http : //www . mpe . mpg . de/erosit a/ 
^ |http : //hea . iki . rssi . ru/SRG/ 
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Even though the evolution of AGN over cosmic time 
is an interesting subject on its own, the remarkable dis- 
coveries of correlations between SMBH masses and host 



galaxy properties (see, e.g., Kormendy & Richstone 1995 
Ferrarese fc Ford||2005[ for reviews), hinting for the under- 



lying AGN-galaxy co-evolution mechanisms, have signifi 
cantly broadened the importance and actuality of those top 
ics in astronomical community. Although controversial in 



terms of observat ional data (e.g., Ferrarese 2002} Kormendy 
& Bender 2011), in currently favored cold dark matter 
(CDM) cosmologies, where the structure grows in a hi- 
erarchical fashion starting from the smallest scales (i.e., 
the bottom- up scenario), it is quite natural that SMBH 
masses A/bh should also (at least up to some level) corre- 
late with the masses of the host dark matter halos Mh. 
Indeed, this is what one quite commonly finds in semi- 
analyt i c AG N-galaxy co-evolution models, e.g. |Fanidakis| 
et al. (2012). Despite the fact that there is a significant 
scatter in the Mbh-M}^ relation, hinting that there are other 
important parameters beyond Mh (and possibly one needs 
a full merger history) that determine Mbh, it is still rea- 
sonable to use Mh as a single proxy for Mqh, this way 
linking Mbh to the quantities we have currently good the- 
oretical predictions available, i.e., halo mass function (MF) 
and clustering bias. 

Future wide field X-ray surveys like eROSITA promise 
to give us very precise AGN clustering bias parameters as a 
function of redshift and luminosity ( Kolodzig et al.|["2013 |. 
Since clustering bias provides direct information on how 
AGN populate dark matter halos, this also helps in deter- 
mining, e.g., the dominant X-ray AGN fuelling mode: 'hot 
halo' vs 'c old gas' accretion. As a n example, semi-analytic 
models of Fanidakis et al. (20121 suggest for the hot halo 
fuelling mode quite strong dependence between the AGN X- 
ray luminosity and host halo mass, with brighter AGN pop- 
ulating preferentially higher mass (and thus more strongly 
biased) halos. For the cold gas accretion, on the other hand, 
which is presumably responsible for most of the quasar ac- 
tivity, the host halo mass and AGN luminosity are only 
weakly dependent, with most of the activity occurring in 
somewhat lower mass halos. 

In this paper we assume that there exists a relation 
between the dark matter halo mass and the AGN X-ray 
luminosity it harbors, and therefore, the X-ray luminosity 
function of AGN can be computed from the mass function 
of dark matter halos. The search for the functional form of 
this relation consistent with the concordance ACDM halo 
mass function and the observed X-ray AGN LF is presented 
in Section [2] It turns out that this way we indeed obtain 
acceptable fit to the X-ray AGN LF and, as a bonus, au- 
tomatically obtain a prediction about X-ray AGN cluster- 
ing properties, which seems to agree with available obser- 
vational data. This and other consequences of our simple 
model is given in Section [3J We bring our conclusions in 
Section g] 

Throughout this paper we assume a flat ACDM cosmol- 
ogy with On, = 0.27, Ob = 0.045, h = 0.70 and erg = 0.80. 



2. Mapping halo MF to X-ray AGN LF 

AGN luminosity can be parameterized via the Eddington 
ratio 
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Fig. 1. Best fitting duty cycle values /dutv('2^i) {upper panel) 
and power law indices T{zi) from Eq. (pi) {lower panel) in 
nine redshift bins, shown with red crosses. In the upper 
panel the power law index is assumed to be in parametric 
form r(2;) = C4 -|- C5Z while the amplitude factors in each 
redshift bins are assumed to be free parameters, i.e., includ- 
ing Mo the fitting problem has 12 free parameters in total. 
With solid black line we show the best fitting three param- 
eter analytic function with the detailed form as specified in 
the legend. In the lower panel /duty (-z) is fixed to the above 
analytic form while the power law indices in each redshift 
bin T{zi) are allowed to vary freely, i.e., this time we have 
13 free parameters. Solid black and gray lines show linear 
and quadratic polynomial approximations to best fitting 
T{zi). The dotted lines and gray shaded regions represent 
best fit values along with la intervals in case all r(zi) and 
/duty(^i) are all allowed to be free (19 free parameters). 



If the Eddington ratio was independent of redshift and 
black hole mass, and the latter was proportional to the 
dark matter halo mass, a simple relation ix oc Afj, would 
obviously hold. In reality, however, the Eddington ratio is 
determined by the environment of the supermassive black 
hole. Although the full picture is likely more complex, one 
can attempt to build a minimalistic model assuming that 
in the population of AGN, on average, the Eddington ratio 
is determined by the dark matter halo mass and redshift, 
/Edd = /Edd(Mh, z) . Similarly, one can assume that the su- 
permassive black hole mass is also on average determined 
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Table 1. Best fitting parameter values for Models I and II. Both constrained and unconstrained cases are given. 
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+ /duty < 1 




Lo = 10« erg/s 
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Cl 


0.47 




0.48 




0.45 


0.43 


C2 


0.83 




0.94 




0.80 


0.90 


C3 


0.11 




0.15 




0.14 


0.21 


C4 


0.84 




0.83 




0.87 


0.87 


C5 


-0.21 




-0.20 




-0.28 


-0.28 


C6 












0.023 


0.023 


Cj= In Mo 


27.3 




27.4 




27.4 


27.5 



Table 2. Global performance of various fitting forms. 



Model 


# param 


# dof 


x' 


expected x 
(lo- CL) 


PLE 

LADE 

LDDE 


7 
8 
9 


72 
71 
70 


96.2 
67.4 
57.8 


72 ±12.0 
71 ±11.9 
70 ±11.8 


Model I 

Model I +/duty < 1 


6 
6 


73 
73 


77.2 
83.8 


73 ±12.1 
73 ±12.1 


Model II 

Model II +/duty < 1 


7 
7 


72 
72 


71.6 
77.9 


72 ±12.0 
72 ±12.0 



by the halo mass and redshift, A/bh = MBuiMy,, z), there- 
fore Lx= J^x(Mh, 0). 

Motivated by these rather trivial and obvious consider- 
ations, we will start by assuming that there is a following 
simple scaling relation between the halo mass and AGN 
X-ray luminosity 



Mu^M, 



Lo 



r(z) 



(2) 



where we take Lq — 10^^ erg/s. Having the above scaling 
relation, we can immediately write for the AGN X-ray LF 



dn 
dZ^ 



{Lx,z) = fduty{z)-r^—{Lx,z) 



dn 
dMh 



[A'ULx,z),z] ,(3) 



where dn/dMh is the concordance ACDM model halo MF 
and /duty is the 'duty cycle', in our case defined as a frac- 
tion of halos that contain AGN in its active state. For the 
halo M F we use the analytic form given by |Sheth fc Tormen 



(T999| and for the X-ra y AGN LF t he ha rd ba nd (2-10 



keV) m easurements from Aird et al. (2010). The Aird et al. 
(2010') LF data points in nine redshift bins arc shown in 



Fig. [3J We assume that the LF measurements in all the lu- 
minosity bins are statistically independent with the errors 
following Gaussian distribution, i.e., the likelihood function 
gets factorized and has a simple analytic form. To sam- 
ple the parametric likelihood function we use Markov chain 
Monte Carlo (MC MC) methods, in particular Metropolis - 
Hastings algorithm ( Metropohs et al.|1953 Hastings|1970 ). 
For finding the maximum likelihood point m parameter 
space we use downhill simplex method as described in Press] 



fduty{zi) obtained from out MCMC calculations. We see 
that T{zi) are quite well fitted by a simple linear form. 

As a next step, we fix T{z) to have a linear redshift 
dependence T{z) — C4 + c^z and allow fdutyizi) to remain 
free parameters. So our free parameters are: Mq, C4, C5, 
f duty {zi) (i = 1 . . . 9), i.e. 12 in total. The best fit values for 
fdutyizi) are shown with red crosses in the upper panel of 
Fig. [ij It turns out that these points can be reasonably well 
approximated with the following analytic form 



/duty (2) 



iz-c2y 



+ cl 



C3; 



(4) 



et al.| ( |1992 ). 



To gam some idea for the acceptable analytical form 
for T{z) and /duty(-z) we initially allow these quantities to 
have free values in all of the redshift bins: T{zi), fduty{zi), 
i = 1 . . .9. Thus, including Mq we have 19 free parameters 
in total. The dotted lines and gray shaded regions in Fig. [l] 
show best fit values along with la regions for T{zi) and 



i.e., a Lorentzian profile multiplied by z and added a con- 
stant offset calj 

To check for the self-consistency of the above analytic 
description we now fix /duty(-z) to have form as given in 
Eq. m , but allow the power law indices in all nine redshift 
bins to be free parameters, i.e., in this case we have 13 free 
parameters: Mq, ci, C2, C3, T{zi) (i = 1 . . .9). The best fit 
values for r(2:j) are plotted as red crosses in the lower panel 
of Fig. TJ where with solid lines we also show the best lin- 
ear and quadratic approximations. Even though the linear 
approximation for T{z) seems to capture the redshift de- 
pendence of the power law index quite well, quadratic form 
with its one extra parameter, as it turns out, is statistically 
well justified. 

Thus, our analytic MF to LF mapping is given by 
Eqs. ([3]), ([2|), (B. In the following we call the model with 
linear (quadratic) T{z) Model I (Model II). Also, we assume 
that there is at most one AGN per dark matter halo which 
is turned on with a probability given by /duty, i-e., in that 
case we have to impose a constraint /duty < 1- However, we 
also considered unconstrained /duty, which can be regarded 
as a case with more than one AGN per halo. It is interest- 
ing that even then /duty does not get much above one (see, 
e.g., the shaded region in the upper panel of Fig. nl). 

In Table [1] we show the best fitting values for aUthe pa- 
rameters of Models I and II, with and without additional 
constraint on /duty The calculated x^ values along with 
expectations are given in Table [2| where we also show the 
results for pure luminosity evolution (PLE) , luminosity and 
density evolution (LADE), and luminosity depe ndent den- 
sity ev olution (LDDE) analytic forms taken from Aird et al 
( 2010 ) . Note that, in order to perform a fair comparison be- 



tween different analytic fitting forms, we have recalculated 



Hiitsi et al. 



(2012h, where we 



Note that in our earlier paper 
also used the MF to LF mapping in the form or Eq. S, we had 
a different analytic form for /duty(z), which turns out to give 
somewhat poorer fit to the observed LF. 
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Fig. 2. Marginalized 2D error ellipses representing la and 2a confidence levels along with ID marginalized probability 
distributions (upper panels in each of the columns) for all the free parameters of unconstrained Model II. The black 
dots and red crosses show the maximum likelihood points for unconstrained and constrained (i.e., /duty < 1) Model II, 
respectively. 



the best fit models using the same fitting machinery as for 
our new models. It can be seen that except PLE all of those 
models provide statistically valid descriptions for the obser- 
vational data. Even though our MF to LF mapping does not 
provide as good fit as LADE or LDDE models, it is quite 
remarkable that this more physical model indeed seems to 
work. As a benefit, once MF to LF mapping is fixed, our 
models make several other predictions, which will be dis- 
cussed in the next section. 

As an example of the resulting parameter uncertain- 
ties we show in Fig. [2] la and 2a marginalized error re- 
gions along with ID probability distribution functions for 
all the seven parameters of unconstrained Model II. The 
black dots and red crosses mark the maximum likelihood 
points for unconstrained and constrained (i.e., /duty < 1) 
cases, respectively. We see that all the model parameters 
seem to be quite well determined. The strongest degener- 
acy is seen between parameters C5 and cg, which might leave 
some doubts about the usefulness of the quadratic term in 
r(z). However, as seen from the results shown in Table [21 
the need for parameter cg is statistically well justified, at 
the confidence level corresponding to '^ 2.5a. 



In Fig. |3] we plot the LFs in nine redshift bins for the 
best fitting PLE, LADE, LDDE models, and for our con- 
strained Mode l II. Original data points with error bars, as 
determined by 
our model slig. 



Aird et al. (2010), are also shown. Note that 
itly but systematically overshoots the high- 



est luminosity data point at large redshifts. 



3. Implications and shortcomings of the model 

In this section we investigate some of the consequences of 
the MF to LF mapping models presented above. 



3.1. Comoving AGN density 

We start by investigating the comoving number density of 
hard band selected X-ray AGN and its dependence on red- 
shift. In the upper panel of Fig. |4] we have plotted comov- 
ing number density for three X-ray luminosity intervals as 
specified in the legend. Here the highest (lowest) group of 
curves correspond to lowest (highest) luminosities. Along 
with constrained Model II results for comparison we have 
also plotted number densities resulting from |Ueda et al.| 
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Fig. 3. Best fitting X-ray AGN LFs in nine redshift bins along with observational data from 

to our constrained Model II fits we also sh ow the best fitting PLE, LADE and LDDE models with the corresponding 

analytic forms taken from Aird et al. (2010). 



( 2003 ) and |Airdet al. ( |2010[ ) LF fits. In case of |Ueda et al. 
( 2003 ) we have used directly their best fit LF para meters. 
The comparison of our results with Aird et al. ( 2010 ) LDDE 
parameterization (which as we explained in the previous 
section was fitted using the same fitting routines as for our 
new models) should give some idea about the possible level 
of uncertainty in n{z), since both of the models give statis- 
tically valid fits to LF data. 

This type of figure is often used to illustrate the anti- 
hierarchical growth of SMBHs (as manifested by the AGN 
activity) where the peak in number density moves towards 
lower redshifts as one considers intri nsically less lumin ous 
AGN. However, at least in the case of Aird et al. (2010) LF 
data, the conclusion about the behavior of the maximum 



in n(z) curve may depend on the particular choice of the 
functional form used to fit the X-ray LF data. Indeed, in 
our LF model, the n{z) maxima for the lowest luminosity 
bins stay almost at the same position, while for the highest 
luminosity bin we see gradual movement towards higher 
redshifts. This does not seem to be entirely incompatible 
wi th the pattern of the original LF data points in Fig. 10 



Aird et al. (2010). Furthermore, in our model the rate of 
change of n{z) has quite strong dependence on luminosity, 
with the lowest luminosity AGN having the fastest rise in 
number density going from high redshifts down to z ^ 1, 
the AGN in the intermediate luminosity bin having only 
very mild increase, the highest luminosity bin show a mild 
drop in the volume density with the redshift. Note that the 
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latter is related to the model behavior at large luminosity 
and redshift, discussed below. 



3.2. A possible shortcoming of tlie model 

A careful examination of Fig.|3]reveals that our model over- 
predicts the numbers of most luminous objects at high red- 
shift (see the three lower panels in Fig. [3]) . This may point 
at the shortcoming of the model. Such a shortcoming could 
be caused by several reasons. One possible cause is the lack 
of proper treatment of the Eddington luminosity limit. The 
Eddington ratio analysis later on in this section suggests 
that this factor may become important in the halo mass 
range log Mh Z 13.5 - 14.5 at the redshift z ^ 2 - 3. It 
may also be an indication of the shortcoming of the un- 
derlying halo mass function, somewhat overpredicting the 
abundance of massive halos at large redshifts. 

Another likely reason may be a failure of the the simple 
power-law scaling between the luminosity and halo mass, 
Eq. (pi), at large halo masses and large redshifts. Indeed, 
due to obvious considerations of time available for the black 
hole growth, the most massive halos at large redshifts may 
harbor somewhat less massive black holes than implied by 
Eq. (pi) . The parameters of the halo mass - luminosity scal- 
ing were determined from the global LF fit which is deter- 
mined by the bulk of the LF data points and is not very 
sensitive to the data at the extremities as the contribution 
of high luminosity points to x^ is relatively small, due to 
their small number and also due to their relatively larger 
uncertainties. 

In the conclusion of this discussion we emphasize that 
the model provides statistically valid global description of 
the X-ray LF data (Table pi). Its (possible) shortcoming at 
the high luminosity end anects only a very small number 
of the most massive halos at large redshifts and does not 
diminish the overall predictive power of the model for the 
bulk of typical AGN, as discussed in the following subsec- 
tions. 



3.3. AGN clustering 

As our model has fixed a scaling relation between the halo 
mass and AGN X-ray luminosity, we have immediately defi- 
nite predictions available for the AGN clustering strengthjj 
In the lower panel of Fig. |4]we compare the clustering bias 
as a function o f redsh i ft with the AGN clust ering data from 
Allevato et al.| (|2011|).|Allevat o et al.| ( |201l"| ) used soft band 
X-ray data from GOSMOS field, which covers 2.13 deg^ 
and contains a total of 780 AGN with available redshifts, 
correspondin g to spectroscopic co mpleteness of ~ 53%. In 
their analysis Allevato et al. (2011 ) applied a magnitude cut 
Jab < 23 and redshift cut z < 4, resulting in 593 objects in 
their final sample, with effective spectroscopic completeness 
~ 65%. 

Since the clustering bias in our model depends on AGN 
luminosity, with more luminous AGN having higher bias 
parameters, in order to make a fair comparison we have to 
adjust the effective flux limit appropriately. Since the selec- 
tion of Allevato et al. (2011) sample was done in the soft 
band, we have to transform their flux limit to the corre- 
sponding limit in the hard band. We do this in a simple 
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Fig. 4. Upper panel: comoving X-ray AGN number den- 
sities as a function of redshift for three 2 — 10 keV band 
luminosity ranges as shown in the legend. The luminos- 
ity is increasing from top to bottom. Along with our con- 
strained Model II we h ave also plotted th e res ults arising 



from t he LF models of Ueda et al. ( 2003 ) and Aird et al 



(|2010|). Lower panel: comparison of bias parameters from 
the constrained Mod el II with the observational data from 
[Allevato et al. ( |2011 ) . The four lines assume different effec- 
tive photon iiiaices (F = 1.4,1.6,1.8,2.0) for the approxi- 
mate transformation between soft and hard X-ray bands. 



way by assuming an effective population-averaged spectral 
index rIf]The results for four different values of F are shown 
in the lower panel of Fig. [4] As in the hard band one starts 
to see more and more sources that are obscured in the soft 
band, the effective spectral index around ~ 1.6 appears to 
be most appropriate. Gonsidering this, we see that except 
for the highest redshift bin, our predictions for the bias 
parameter seem to agree reasonably well with the obser- 
vational measurements. The highest redshift bin is repre- 
sented by one point only and it remains to be seen, how 
signiflcant the deviation is. Also, one has to keep in mind 
that the measurement of the linear bias parameter, and es- 
pecially the proper estimation of the errors due to cosmic 



^ In de tail, for the halo c lustering bias we use the analytical 
model oflShethet al.|(|2001|. 



Not to be confused with the power law index r{z) appearing 
in Eq. ^. 
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variance (e.g., accidental superclusters in small fields could 
significantly influence the results) from the fields as small 
as ~ 2 deg^ may be quite a cumbersome task. 

3.4. AGN duty cycle 

In both constrained and unconstrained models the AGN 
duty cycle (Eq. ([3|) peaks at the redshift z ~ 1 (Fig. [1]). 
Although unconstrained models (i.e, when the duty cycle 
was allowed to take any positive value) result in somewhat 
better description of the observed X-ray LF data, it is in- 
teresting that the peak value of the duty cycle in these 
models is nevertheless close to unity, with the upper limit 
of /duty ^ 1-5 (Icr confidence). This implies that on average 
there is at most ~one AGN per dark matter halo. This also 
implies that the volume densities of dark matter halos and 
AGN at the redshift 2: ~ 1 match, i.e. there is enough dark 
matter halos to explain all AGN. In both constrained and 
unconstrained models the duty cycle drops to /duty ~ 0.2 
at the redshift z = (Fig. [l]). 

3.5. Mh — Lx scaling relation 

As stated above, the crucial component in our model is the 
assumption about existence of a scaling relation between 



the halo mass and AGN X-ray luminosity. [Koutoulidis et al. 



(2012) have recently attempted to infer this type of scal- 
ing by using AGN clustering data from several deep X-ray 
fields: Chandra Deep Field South and North, the AEGIS, 
the extended Chandra Deep Field South, and the COSMOS 
field. In Fig. [5] we show their measurements with error bars 
along with our constrained Model II predictions for red- 
shifts 0.5, 1.0 and 1.5, i.e. approximately covering the red- 
shift range of AGN used in their analysis. The dotted line 
with surrounding gray band displays the results from the 
semianalytic galaxy- AGN co-evolution model of [Fanidakis 
et al.| ( 2012[ ) (as presented in Fig. 9 of Koutoulidis et al. 



( |2012[ )j, showing their 'all AGN' case, i.e., the sum of 'hot 
halo' and 'starburst' modes. Although the uncertainties are 
rather large, we see that our model is able to capture the 
major trend where the more luminous AGN populate pref- 
erentially more massive halos. The quantitative agreement 
between ou r M^ — Lx relation and results of |Koutoulidis| 



et al. ( 2012 ) is truly remarkable as they are based on totally 
independent arguments. 

The long-dashed line in Fig. [5] shows the M^^-Lx scal- 
ing one obtains when converting the Mji-Mbh relation of 
Bandara et al. (20091 by assuming a constant Eddington 
ratio /Edd = Lx/LE ddJMBu) of 3 x 10""*. In reality, as 



Bandara et al. ( 2009 1 Mbh-M^ scaling relation used strong 
lensing masses derived from the SLAGS lens sample with 
mean redshift z ~ 0.2, we have adjusted the amplitude 
of the scaling relation to correspond to w hat one might 
find at redshifts z ^ 1, using results from Croton (2009) 
semi-analytic models. In Bandara et al. ( 200^[ ) theauthors 
obtain A/bh oc M^ with 7 = 1.55 ± 0.31, which agrees 
quite we ll with analy tic and s e mi-an alytic models of |Wyithe| 



7 



573 



& Loeb| (|2003]) and 
1.67 anc 



Croton (2009), which provide values 



and Wyithe & Loeb 

scaling relation increases by factors of ' 



1.39, respectively. In |Croton| ( |2009[ ) 
(2003) models the amplitude of the 
2 and ~ 4, respec- 
tively, while going from z = to z = 1. However, the power 
law index 7 is independent of redshift in both models. 
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Fig. 5. Comparison of the Mi^-Lx scaling relation with 



the observational data from Koutoulidis et al. 



(2012). To 



capture the redshift ra nge of majority of AGN used in 



Koutoulidis et al. ( 2012 ) analysis we have plotted the model 
scaling relations for three redshifts: z = 0.5,1.0,1.5. The 
dotted line with surrounding gray band displays the results 
from the semianalytic galaxy- AGN co- evolution model of 
Fanidakis et al. (2012), as presented in Koutouhdis et al. 
]2012), showing their all AGN' case, i.e., the sum of 'hot 
halo' and 'starburst' modes. The long-dashed line shows the 
M\i-Lx scaling one obtains o nce converting the amp litude- 
adjusted Mh-AfeH relation of Bandara et al. ( 2009 1 by us- 
ing the fixed Eddington ratio /Edd = -'^x/^Edd(AfBH) = 
3 X lO-''. 



3.6. AGN Black hole masses, Eddington ratios and the 
growth ofSMBHs 



By assuming the A/eH-A'/h scaling of|Bandara et al.' ( |2009 ) 
with amplitude adjusted according to jCroton ( 2009) semi- 
analytic model, we can convert our LxAfh relation to /Edd- 
Afh scaling relation. This is detailed in Fig. [6] where we have 
plotted various scalings for five different redshifts. The up- 
per panel gives our results for Lx-Afh with la uncertain- 
ties as determined in Section [2] for constrained Model II. 
The m iddle panel shows amplitude-adjusted [Bandara et aL] 
(2009) AfBH-A^h relation with Icr uncertainties. In the lower 
panel we show the Eddington ratios as determined by com- 
bining the results from the two uppermost panels. The X- 
ray luminosities and Eddington ratios refer to the 2-10 keV 
band and no bolometric corrections were applied. For a typ- 
ical type I AGN spectrum, the bolometric correction factor 
should be in the ~ 5 — 10 range. 

As we can see, the Eddington ratios for high redshift 
(z > 2) AGN are predicted to be steeply rising functions of 
halo mass, while at redshifts z ^ 1, i.e. close to the peak 
in X-ray AGN number density, the dependence of /Edd on 
Afh flattens out, and at lower redshifts the trend is mildly 
reversed, i.e., the AGN in less massive halos have some- 
what higher Eddington ratios. Note that for z ^ 3, the 
2-10 keV band /Edd reaches values of ^ 0.1 for halo masses 
~ 2 X lO^'^ Mq, i.e. the bolometric Eddington ratio will ap- 
proach ~ 1 at this mass. At larger masses, our Afh — Lx 
scaling relation would imply /Edd ^ 1- Obviously, this is a 
consequence of the fact that our Afh — Lx scaling relation 
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Fig. 6. Upper panel: average Lx-M^ relations (Eq. ([2|) 
along with la uncertainties for five different redshifts 
as given in the legend. Middle panel: amplitude-adjusted 
MBH-Mh sca ling relation with la uncertainty from |Bandara| 
et al. ( 2009 1 . Lower panel: Eddington ratios in 2 — 10 keV 
band for hve redshifts arising from Lx-My^ and M^h-Mi^ 
relations given in the upper two panels. 



does not impose the Eddington luminosity limit. This may 
be considered as one of the shortcomings of the model, as 
it was discussed in the section [T2l One has to note how- 
ever, that this does not affect overall performance of the 
model, in the major part of the parameter space. Indeed, 
so massive halos at redshifts z ^ 3 are extremely rare, cor- 
responding to ~ 4 — 5tT peaks in the initial fluctuation field. 
Bulk of the objects are shining at much smaller Eddington 
ratios, <C 1. 

Having specified Mbh-M^ scaling relation we can also 
calculate how SMBHs with different masses grow through- 
out the cosmic time. Taking into account the duty cycle, 
the population average e-folding time of the SMBH growth 



T = 



Mbh 



V 



Mbhc^ 



M /duty /duty /bol Lx 

where rj is the accretion efficiency, defined as 



L. 



fboiVMc' 



(5) 



(6) 



p 




Fig. 7. The hard band X-ray luminosities {upper panel) and 
SMBH e-folding times r = ^^^^ in units of the Hubble time 

(lower panel) for various black hole masses as shown in the 
legend. For the low mass SMBHs at higher redshifts, where 
they are too faint to be constrained by the observational 
data, we have shown with dotted lines the extrapolations 
of the model curves and have omitted the corresponding la 
uncertainty bands. 

Lx refers to the X-ray luminosity in the 2-10 keV band used 
throughout this paper and /boi ~ 5 — 10 is the bolometric 
correction factor for the 2-10 keV band. The e-folding time 
can be further expressed though the Eddington ratio: 



77 Mbhc^ 



/duty /bol /Edd L 



= 0.45 



'/ 



Edd 



/duty /bol J] 



Edd 



Gyrs (7) 



In the lower panel of Fig. [7] we plot the e-folding time in 
units of the Hubble time, i.e. tH{z), vs redshift for several 
values of the black hole mass. Shown in the upper panel is 
the evolution of the population average 2 — 10 keV band 
X-ray luminosity with redshift. AGN harboring black holes 
with masses less than a ^fewx 10^ AIq are too faint at large 
redshifts a nd fall outside the luminosity range of the X-ray 
LF data of Aird et al. (20101 used in this study. Therefore 
they are not directly constrained by the X-ray data. At 
these redshifts we only show with dotted lines the extrapo- 
lations of the central model curves and omit the la uncer- 
tainty regions. 

We see that the X-ray luminosity of more massive 
SMBHs declines faster with the redshift than for less mas- 
sive ones. For ^ 10^ Mq black holes Lx stays almost con- 
stant while for even lower mass SMBHs the evolutionary 
trend is slightly reversed, i.e., the objects are getting some- 
what brighter as one moves towards lower redshifts. As one 
might have expected from the results plotted in the bottom 
panel of Fig. |6j for the typical X-ray selected AGN the cor- 
responding SMBH e-folding times are large in comparison 
to the Hubble time. With the exception of most massive 
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SMBHs at large redshifts, the typical X-ray selected AGN 
are not in the rapid mass growth regime in the z ^ — 3 
redshift range. 

Note that the crossing of the model curves in the lower 
panel of Fig. [7] at the redshift of z ~ 0.85 is artificial and 
does not seem to have any physical meaning. It is driven by 
the equality of the power law indices of the Mbh-My^ and 
M\i — Lx scaling relations, which leads to the tH being 
independent of Mbh at that particular z. 

Thus, wit h the particular form of the M^i^-Mi^ scaling 
relation fr om Bandara et al. (2009) with amplitude evolu- 
tion from Croton (2009), our model implies relatively low 
Eddington ratios, jEdd ~ 10^"* — 10^"^. In the redshift range 
z ~ — 3, for a typical X-ray selected AGN accretion pro- 
ceeds in the 'hot halo mode' characterized by relatively low 
mass accretion rate and large e-folding time for SMBH mass 
growth. Obviously, this conclusion critically depends on the 
assumed MBn-Mh scaling relation. Indeed, as X-ray lumi- 
nosity for the given halo mass is directly determined by 
our model from the Lx — M^ mapping, the Eddington ra- 
tio scales inversely with the black hole mass, while e-folding 
time is linearly proportional to the black hole mass. 



4. Conclusions 

Even though dark matter halos are highly nonlinear ob- 
jects, with the aid of cosmological N-body simulations along 
with heuristic analytic models, one has obtained a very 
good knowledge about their redshift dependent MF and 
clustering bias. Since in CDM type cosmologies, where the 
structure forms according to the 'bottom-up' scenario, it is 
quite natural to assume that the buildup of more massive 
object has also been accompanied by the formation of more 
massive SMBH in its center. This is also manifested on av- 
erage by a more luminous AGN. Despite the fact that there 
are certainly several other key parameters beyond the halo 
mass, which influence the central SMBH mass, Afh is surely 
expected to be one of the main contributors. 

In this paper, we assumed that a scaling relation be- 
tween X-ray AGN luminosity and its host halo mass does 
exist. Assuming further that the dark matter halo mass 
function is described by the concordance ACDM model, we 
can predict X-ray AGN luminosity function. Comparing the 
latter with the redshift dependent AGN X-ray luminosity 
function known from observations we can determine the 
shape and parameters of the Mji — Lx scaling relation. 

Our main conclusions are the following: 

— We have shown that a simple scaling relation in the form 
of Eq. ([2|, where T{z) is a linear or quadratic polyno- 
mial, allows one to obtain X-ray AGN LFs consistent 
with observational data. The best fitting model param- 
eters along with x^ values are summarized in Tables [I] 
andO 

— Since the Mi^-Lx scaling establishes a link to the halo 
MF, our model is certainly more predictive compared 
to the usual specifically tailored multi-parameter fitting 
forms to the X-ray LF data. In particular, our model 
predicts the redshift and luminosity dependence of the 
X-ray AGN clustering bias, which is in reasonable agree- 
ment with observational data. 

— Our Afh-Lx scaling rel ation is in good agr e ement with 
the data presented in Koutoulidis et al. (2012), and 
with semianalytic galaxy-AGN co-evolution models of 



Fanidakis et al. (2012). Comparison with semianalyti- 
cal models suggests that for X-ray AGN the dominant 
accretion mode is the 'hot halo mode', in contrast to 
the 'starburst mode' which is compatible of being dom- 
inant accretion mode for somewhat less clustered but 
more energetic quasars. 

Our model also includes AGN duty cycle. It peaks at 
the redshift z ~ 1 and its value at the peak is consistent 
with unity, implying that there is at most one AGN per 
dark matter halo. 

We further co mbined our Mh - Lx re lation with Mh-MeH 
relation from 'Bandara et al.'('2009), along with its red- 
shift evolution taken from Croton ( ^2009) , and obtained 
Eddington ratios /sdd = ^x/^Edd ^'^^Hlr^ _ ;^q-3 (•2_io 
keV band, no bolometric correction applied). This result 
presents another evidence that at the redshifts below 
z ^ 3, SMBHs in X-ray selected AGN typically grow 
through 'hot halo' accretion mode characterized by low 
mass accretion rate and large e-folding time. It is also 
consistent with the fact that according to their cluster- 
ing strengths. X-ray selected AGN populate group-sized 
dark matter halos and thus are typically embedded in 
hot gaseous halos. 



Appendix A: Analytic fits for concordance ACDM 
halo MF and bias parameters 

To facilitate a fast check of our results for those who do 
not have numerical routines for the halo MF and cluster- 
ing bias available, we provide here analytic fits for those 
quantities, valid for the mass and redshift ranges relevant 
for this study: A4 > 10^^ Mq, z <3~A. 
Halo MF can be approximated as 



dn 



dMh 

where 

A{z) 
]VL{z) 



(Afh, z)~A(z) 



AL{z) 



a{z) 



exp 



A4 
A'L{z) 



P(z) 



exp (-0.1057 -z^ 



2.677 • z - 43.89) , 

exp (0.08305 • z^ - 1.709 • z + 33.17) , 

a{z) = 0.007189 • z^ - 0.1692 • z - 1.919 , 

/3(z) = 0.001109 -z^- 0.04119 •z + 0.7188. 

Similarly, for the halo clustering bias, as a function of halo 
mass and redshift, one approximately obtains 

&(Mh, z) ~ ci(Afh) • z2 + C2(A'fh) • z + C3(Afh) , 

where 

ci(A4) ^ exp [0.017031n(Afh)2- 0.6629 ln(A4) 

C2(A4) = exp [0.014951n(Mh)2- 0.5196 In(Mh) 

C3(A4) = exp[0.028661n(A4)2-1.5341n(A/4)4 

Under the mass and redshift constraints as given above, 
those analytic f orms provide fits to the |Sheth fc Tormen| 
( |1999[ ) MF and jSheth et al.| ( |2001|) bias parameter with 
accuracies ^ 10% and better than 5%, respectively. NOTE 
that in the above formulae we have assumed that MF is 
measured in units of /i'^Mpc~'^/ft.~^Af0 and halo mass in 
/i~^Mq, while in the main part of our paper we fixed the 
reduced Hubble parameter to h = 0.7. 
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